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Qh! Abstract 
^ ■ 

• • . We complete the investigation of loop corrections to hard spectator-scattering in 

I exclusive B meson to light meson transitions by computing the short-distance 

^ ' coefficient (jet-function) from the hard-collinear scale. Adding together the two 

. coefficients from matching QCD — > SCETj — > SCETn, we investigate the size of 

loop effects on the ratios of heavy-to-light meson form factors at large recoil. We 
find the corrections from the hard and hard-collinear scales to be of approximately 
the same size, and significant, but the perturbative expansions appear to be well- 
behaved. Our calculation provides a non-trivial verification of the factorization 
arguments. We observe considerable differences between the predictions based on 
factorization in the heavy-quark limit and current QCD sum rule calculations of 
the form factors. We also include the hard-collinear correction in the B — > vrvr tree 
amplitudes, and find an enhancement of the colour-suppressed amplitude relative 
to the colour-allowed amplitude. 



1 Introduction 



The matrix elements of flavour-changing currents q Tib are important strong interaction 
parameters in low-energy weak-interaction processes. The strong interaction dynamics 
of semi-leptonic B decays is encoded in these form factors. They are also inputs to the 
factorization formulae for hadronic two-body B decays P and radiative decays [2]. A 
better understanding of such quantities improves the accuracy of the extraction of the 
CKM matrix parameters from experimental data, and of searches for new phenomena in 
flavour- changing processes. Thus efforts are being made to compute the form factors with 
different methods including QCD lattice simulations [HI, light-cone QCD sum rules jlj 
and quark models 0. 

It is also interesting to investigate these form factors in the heavy-quark expansion. 
It is well-known that all B ^ D^*^ form factors reduce to a single (Isgur-Wise) function 
[H] up to calculable short-distance corrections at leading order in this expansion. In this 
paper we consider transitions of B mesons to light mesons in the large-recoil regime, 
where the light meson momentum is parametrically of order of the heavy-quark mass. 
In this regime a similar simplification applies to heavy-to-light form factors [7j: the 
three (seven) independent B pseudoscalar (vector) meson form factors reduce to 
one (two) function(s) up to corrections that can be calculated in the hard-scattering 
formalism at leading order in the heavy-quark expansion [S]. The different form factors 
can therefore be related in a systematic way. The factorization formula that summarizes 
these statements reads jH] 

Ff'-^'iE) = C,{E)UE)+ / — / dvT,{EM^,v)ct>B^{uj)ct>M{v) (1) 

JO UJ Jo 

with E the energy of the light meson M, ^a(-E') the single non-perturbative form factor 
(one of the two form factors when M is a vector meson), and (px the light-cone distribu- 
tion amplitudes of the B meson and the light meson. The short-distance coefficients Ci 
and the hard-scattering kernel Tj can be calculated in perturbation theory. The heavy- 
to-light form factors are more complicated than both, the B — > D^*^ form factors and 
light-light meson transition form factors at large momentum transfer. Contrary to the 
case of i? ^ D^*\ a spectator-scattering correction, the second term on the right hand 
side of (Q), appears. On the other hand, the form factor cannot be expressed in terms of 
a convolution of light-cone distribution amplitudes alone, because the corresponding con- 
volution integrals are dominated by endpoint singularities jU] . In ([Q) these contributions 
are factored into the the function ^a{E).^ The factorization formula ([Q) has been shown 
to be valid to all orders in perturbation theory ^21 (see also ^3 EI) in the framework 
of soft-coUinear effective theory (SCET) [T31 Uni Ej . In particular, since the two rele- 
vant short-distance scales mi, and (m^A)^/^ (A is the characteristic scale of QCD) can be 

^The statement that the endpoint contributions are not calculable is challenged in the PQCD ap- 
proach |10| . which assumes that Sudakov resummation renders them perturbative. This point is critically 
examined in . We also note here that our notation (E) does not show the dependence of the form 
factor on the nature of the meson M. 
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separated in SCET, the short- distance coefficients Tj pertaining to spectator-scattering 
are represented as convolutions C^^^^ -k J with the two factors associated with the two 
different scales. 

In the limit that not only power corrections in A/m^ but also radiative corrections 
in the strong coupling are neglected, the second term on the right hand side of (P) 
is absent, and parameter-free relations between ratios of form factors follow [7]. The as 
contributions to (H} have been computed in [H], and the spectator-scattering term Tj has 
been found to dominate the correction. This motivates an investigation of the subsequent 
term in the perturbative expansion of Tj. Since the leading term is due to a tree 
diagram with gluon exchange between the current quarks and the spectator anti-quark, 
this amounts to the computation of the 1-loop correction to spectator-scattering. Since 
Tj = Ci^^^-k Ja, the calculation splits into two parts. In a previous paper [TH] (see also 
|19j ) we reported the ffist part of the calculation which consisted of the 1-loop correction 
to the coefficients originating from the hard scale m^. In this paper we complete 

the calculation with the 1-loop computation of the "jet-functions" Ja originating from 
the "hard-collinear" scale (mfeA)^/^. The jet-functions have also been computed by Hill 
et al. fniEDj- Nevertheless, an independent calculation is useful, since the computation 
is quite involved and a comparison showed that the result of [20] was originally not given 
in a scheme consistent with the MS definition of the light-meson light-cone distribution 
amplitude (see the discussion in ^nillOl)- Furthermore, the numerical impact of these 
calculations on the relation between form factors and other observables in B decays has 
not yet been discussed in any detail in the literature. 

The organization of the calculation of the short- distance coefficients Ci and Tj follows 
closely the derivation of the factorization formula in J2|- In a first step, the effects from 
the hard scale are computed and QCD is matched to an intermediate effective theory, 
called SCETj. In SCETi the term and the hard-scattering term are naturally defined 
by the matrix elements of two distinct operator structures, the so-called yl-type and 
-B-type operators. At this step, the form factors can be represented as 

F^-'\E) = Q{E)UE) + / dTCr\E,T)Ea{T,E). (2) 

The point to note here is that the three (seven) form factors oi a. B P {B ^ V) 
transition can be expressed in terms of one (two) form factor (s) C,a{E) and one (two) 
non-local form factor(s) 'Ea{T,E). A number of relations between form factors emerge 
already at this stage. In Section |21 we define the SCETi operator basis, and express the 
QCD heavy-to-light form factors in terms of the SCETi hadronic matrix elements, which 
leads to (j21). All the required short-distance coefficients of the SCETi operators can be 
inferred from ^Hl ^] • 

Eq. ^ is useful only to a limited extent, because it introduces the form factors 
5a(r, £^), which depend on two variables. However, it has been shown that, contrary 
to the ^a{E), the Sa(r, i?) can be factorized further into a convolution of light-cone 
distribution amplitudes with a hard-scattering kernel (jet-function) jT2]. This amounts 
to performing a second matching to SCETn, in which the effects at the hard-collinear 
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scale (mfeA)^/^ are computed. This is done in Section |21 Here we discuss in detail the 1- 
loop calculation and renormalization of the jet-functions Ja that follow from representing 
the SCETi matrix element of the B-type operators in the form 

Ea{T,E) = - / du dv Ja{T]V,U}) fB4>B+{'^)fM4>M{v). (3) 

4 JO Jo 

Combining this with (0) we obtain the spectator-scattering term in (^. The calcula- 
tion is done in dimensional regularization which requires dealing with evanescent Dirac 
structures specific to d dimensions. As will be discussed, a subtlety arises due to the fact 
that the factorization properties of SCETn require a specific choice of reduction scheme. 
Together with Ja we also determine the anomalous dimensions of the B-type operators 
confirming the results of |2Uj . 

The detailed numerical analysis of the corrections from the two matching steps is 
contained in Sections |^ and In addition to the next-to-leading order correction we also 
include the summation of formally large logarithms from the ratio of the hard and hard- 
coUinear scale by deriving a renormalization group improved expression for the coefficient 
functions Cl . From the size of the 1-loop correction we conclude that the perturbative 
calculation of spectator-scattering is under reasonable control despite the comparatively 
low scale of order (m;,A)^/^ ~ 1.5 GeV. The combined hard and hard-coUinear 1-loop 
correction is about (50 — 70)% depending on the observable. This is also of interest in 
the context of QCD factorization calculations of hadronic B decays, since the same jet- 
function enters the spectator-scattering contributions to two-body decays [25 . Sectional 
is devoted to a discussion of the symmetry-breaking effects on the form factor ratios 
and a comparison of these ratios to QCD sum rule calculations. We then consider the 
tensor-to-(axial-)vector form factor ratios that appear in electromagnetic and electroweak 
penguin decays, and the numerical impact of our jet-function calculation on hadronic 
decays to two pions. Here we find that the new contribution increases the ratio of the 
colour-suppressed to the colour-allowed tree decay amplitude, which leads to a better 
description of the branching fraction data. We conclude in Section El 

In Appendix^ we summarize the short-distance coefficients C\^^\e,t) relevant to 
(121). Some of the convolution integrals of the jet-functions and the coefficients C\^^^ {E, r) 
needed for the numerical analysis of the spectator-scattering term are collected in Ap- 
pendix ini 

2 Heavy-to-light form factors in SCETj 

Our first task is to express the QCD form factors in terms of matrix elements of SCETi 
currents and the corresponding short- distance coefficients. We use the position-space 
SCET formalism and the notation of ^21 QBl to which we refer for further details. 
The "collinear" fields ^ and Ac that appear in this section describe both, hard-coUinear 
(virtuality rribA) and collinear (virtuality A^) modes. The reference vectors v, Uzp are 
defined such that v"^ = 1, = = 0, n_n^ = 2. Except for Section ITT] we adopt a 
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frame of reference where n^v = 1 and v = (?7,_ + n+)/2. In scalar products of n_, n+ 
with other vectors we omit the scalar-product "dot" . 



2.1 Operator basis 

The relevant terms in the SCETj expansion of a heavy-to-light current ipTiQ read ^3 

dllHl 



(vr,Q)(o) = j dsY.c\f'\s)of'\s-Q) 
j 

j 



where 



Of^^\si, 32-, x) = — {^W,)sAW^,iD^,^W,)s,^K, (5) 



1 

rrib 

and Si = Sirrih/n^v. Since the coUinear fields ^ and Ac describe modes of different virtu- 
ality, no simple A/mfo-scaling rules apply to these fields. The power-counting argument 
that shows that the three types of operators contribute to the form factors at leading 
power in the heavy-quark expansion has been given in J2] • The main difference between 
the two types of operators is their dependence on position arguments. The B-type op- 
erators are tri-local, and for this reason are sometimes also referred to as "three-body" 
operators. The 1-loop corrections to the coefficient functions of the A-type currents have 
been calculated in [121 EE] , to those of the B-type currents in [THl CH] • 

The basis (0) is motivated by the simple expressions of the tree-level matching coef- 
ficients in this basis ^Hl- However, the analysis of ^2] shows that O]^^'' and Oj^^^ are 
operators relevant at leading power in the 1/mb-expansion only because of the trans- 
verse coUinear gluon field in the covariant derivative D±c- It is therefore advantageous 
to perform a basis redefinition such that the transverse coUinear gluon field appears only 
in Oj^^\ This can be done by replacing Oj^^^ by 

{^Wc)stt^^{tn_vn+t)"%K, (6) 

which is the choice that has been adopted in 1^ . The redefinition involves the 
identity 
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= {^w,)s ' T'-K - {^w,)s{w^iD^,,w,), ^ r'^K (7) 

The second term modifies C\fi^\si,S2) in the new basis by an amount proportional to 

C'lf^\s)? In the new basis only the AO- and B-type operators contribute to the form 
factors at leading power in the 1/mb-expansion. Our new basis of operators for a given 
Dirac structure Vi is: 



Scalar current J = ipQ: 



\ tn, a / 



j(Bi) ^ p^cW,]K (8) 

Vector current J^^ = ip'~^^Q: 

\ m+d ^ y n-v n_v d 
^ mi, n_v 

jiBiH ^ ^{^W,)^,, [W}i ip^mK (9) 
Tensor current J^^, = ipia^uQ- 



^In momentum space the generic modification is c[^^\e,t) — cj^f^\E,T) +mb/ {2E) C^f^\E), see 
also Appendix IXI 
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/ / i^f It — (J 

4f = —m^h[,.i^.M^ P^cW^] K + Jif )^ (10) 

Here a[^6y] = a^fe^ — a^6^. The operator J^^^^"^ vanishes in four dimensions, but 
must be kept since we regularize dimensionally. 

Here we dropped the position indices si,2 which should be clear from (0). We also 
dropped the operators involving explicit factors of position x^, which come from the 
multipole expansion (see JHl); since we can always work with the QCD currents at a; = 0. 
The choice of the AO operators is identical to that of [201; t)ut the basis of B-operators 
is slightly different. As in [201 we combined the Al-operators with the AO-operators 
using that their coefficient functions are related |TS1I22I- Since the Al-operators without 
the transverse hard-coUinear gluon field do not contribute to the form factors at leading 
power ^2]; these extra terms to the J'^^^-' will not be considered in the following. The 
SCET representation of the QCD current Jy(0) is then 

^x(O) = E c^f^' * Jx'^' + E c*?' ^ Jx'^' + ■ ■ ■ > (11) 

i k 

which defines the coefficient functions for the scalar (X = S*), pseudoscalar (P), vector 
iy), axial- vector [A) and tensor (T) currents. The star-product of coefficient function 
and operator in position space is a convolution over the arguments Si as in Pj). 

The basis of the pseudoscalar (axial-vector) operators can be inferred from the scalar 
(vector) basis by the replacement {^W^ — > (^PVc)75. In a renormalization with anti- 
commuting 75 (as adopted in ^H] and in this paper), the short- distance coefficients 
q{M)i ^ (j{Bi)k ^j^^ scalar and pseudoscalar current are then equal, as are those of the 
vector and the axial- vector current J^s = '?^757^Q (note the order of 75 and 7^). In 
Appendix ^ we give the transformation of the momentum-space coefficient functions 
calculated in ^H] to the new basis. For the remainder of the paper we adopt the frame 
where n_w = 1 and v = (n_ -|- n+)/2. 



2.2 Definition of the QCD form factors 

The matrix elements of the QCD currents are decomposed into Lorentz-invariant form 
factors. Following the conventions of [H] the independent form factors are 



{P{p')\qYb\B{p)) = Uiq' 



— m 



+ /o(g 
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{Pip')\qa'^''qub\Bip)) 



rriB + mp 



q^{p^ + p'^) - {mg - mp)q 



2 



for pseudoscalar mesons, and 



{V{p\e*)\q^^^h\B{p)) 



2iV{q^) 

TUB + my 



e'"''''e:p'p., 



{V{p',e*)\qr75b\B{p)) = 2myA,{q'f—^q'^ + imn + my)A,iq') 



-A2{q' 



q 

e* ■ q 



'p'' + p'^ - 

g2 



TTlB + niv 

{V{p',e*)\qta^''q,b\B{p)) = 2zTi(g2)e/^-P'^e»., 
{V{p',e*)\qta^''^5q.b\B{p)) = T2(g2)[(m| - m'y)e*'^ - {e* ■ qW + p'^ 



(12) 



q^ 



+Uq'){e*-q) 



q' 



m-jj - my 



-{p^+p'^) 



(13) 



for vector mesons. We define q = p — p' and use the convention e°^^^ = — 1. From now 
on we neglect terms quadratic in the light meson masses mpy, but keep linear terms. 
In this approximation we can put p' = En^ with E = n^p'/2 = {m\ — q^)/(2mB), and 
e* -71- = 0. 



2.3 Definition of the SCETi form factors 

Taking into account the quantum numbers of the mesons, it is straightforward to relate 
the matrix elements of the operators defined in (jHJ to ()10|) and the corresponding pseu- 
doscalar and axial-vector operators to a few non-vanishing SCETi matrix elements. We 
first note that one of the Sj-integrations in ^ can be done explicitly using collinear mo- 
mentum conservation ^21- This allows us to focus on matrix elements of AO-operators 
with s = and of B-type operators with si = 0. To see this for the case of the B-type 
operators, we represent the position-space coefficient functions in terms of 

where the arguments Xi of the momentum-space coefficient functions correspond to 
the momentum fractions Xj = n^pyrrif, of the collinear building blocks {^Wc)si and 
(W^iD'^^Wc) S2 of the current operator. Then with (jH) and (0) we obtain 

(M(p')l / dhdhCg'\h,h)Of'^''is,,s,-0)\B^) = ^Jd^Cg'^{^^^) 

X {2E) I ^ e-^'^^^ (M(p')|(eiVc)(0)(iyi^Z?::,l^.)(rn+)r;./i.(0)|5.) (15) 
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with f = 1 — r. Abusing notation we will write the coefficient functions simply as 
C\fil\E,T) in the following. Next, the J'^'^^^-operators (except for the tensor operators 
with two transverse indices, which we do not need in the following) can be written as 
(Fourier transforms of) linear combinations of 

Mr) = 2E I {lw:)mwliD^^,W:){rn+)TuhM, (16) 

J Zti 

where for k = 1, 2, 3 the Dirac matrix can take one of the three expressions 

Tfe = {(75)7^S (75)7.x7"S (75)7'^"7.x}, (17) 

and f = rrrih. Here the 75 in brackets means that this factor may be added. This notation 
is convenient, because many of the results below do not depend on the extra factor of 
75. In the following definitions we leave out the position argument of a field, when it is 
X = 0. Eq. (fTH|l suggests defining the B-type form factors as the matrix elements of the 
operators Jk{T). 

We therefore define the two leading-power SCETi form factors for pseudoscalar 
mesons through 

{P{p'mW,)K\B,) = 2E^p{E), 

{P{p'mW,){Wliip,^W,){rn+)K\B,) = 2m^E j dr e'^''^'' Ep{t, E). (18) 

Here \By) denotes the B meson state in the static limit (see the Lagrangian (j32|l below) 
normalized to 2m b (rather than 1 as is conventional in heavy quark effective theory). 
The second definition is such that 

Ep{T,E) = i2m,E)-' {P{p')\Mr)\B,) (19) 

(no 75 in i7i(r)). Similarly, for vector mesons 

{V{p'mW,)j,K\B,) = -2Ee*-v^\\{E), 

{V{p'mWcHlf.MB,) = -2E {e^-e-vn.^)UE), 

{Vip'mW,h5{W^iip,^W,){rn+)K\B,) 

= -2mbEe*-v J dr e'^^^' E\i{t, E), 

{Vip'mWch5l>.±{W^c^iPc±W,){rn+)K\B,) 

= -2m,E [el-e*-vn^^) j dr e'^''^'- E^{t, E), 

(^(p')l(eX)75(W^i^^c±W^c)(rn+)7^^/i,|5,) 

= -2m,E{e;-e*-vn_^) J dr e'^''^'- E^{t, E). (20) 
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The tensor operators J^f^^^, J^^^^'^ must have vanishing matrix elements between a 
pseudoscalar B meson and a pseudoscalar or vector meson, so the set of B-type operators 
J7i_3(r) (now including 75) is complete. We shall find in SectionElthat the matrix element 
that defines H_|_(r, £') (corresponding to J^{t)) vanishes at leading order in the 
expansion, hence we set Sx(r, ii^) = in the remainder of this section. The dependence 
on the polarization vector e* shows that the form factors with subscript a = || (a = ±) 
refer to longitudinally (transversely) polarized vector mesons. 

2.4 Form factor expressions 

Taking the matrix element of (fTT|) . using (fTH|l . and inserting the definitions of the QCD 
and the SCETj form factors we derive expressions for the QCD form factors, in which 
the effects at the scale are explicitly factorized. For the pseudoscalar meson form 
factors, we find 

U{E)=C\f\E)ip{E) + J drCf^'\E,r)Ep{T,E), 
fo{E) = djC'\E) UE) + J drClf'\E,T)Ep{T,E), 



ME) = C\f\E)UE) + J drC\f\E,r)Ep{T,E). (21) 



iriB + mp 

Similarly, for the form factors of vector mesons 

"^l" V{E) = c'y^'\E) U{E) + I dTdy'''\E,T)E^iT,E). 
niB + my J 

^ME) = djy{E)^,iE) + I drd,^'\E,r)E,{r,E), 



2E 



A,{E) = dy'''\E)UE) + J dTC{f'\E,T)E^{T,E) 



^^^^^AiE)-^^^^^ME) 
IE Trip 

= Ci^'\E) ^\\{E) + J dTC\f\E,r)E\\{r,E), 
n{E) = d^y{E)UE) + I dTC^^'\E,T)E^{T,E), 
^T,{E) = Ci^;\E)UE) + J dTC!f^'\E,T)E^{T,E), 
^ (E) - T3 (E) = Cjf ) (E) (E) + JdT C}f ) (E, r) H|| (r, E) . (22) 

We recall that E denotes the energy of the light meson. The coefficient functions C^^^"^ 
and dp^^^ are defined as linear combinations of momentum-space coefficients functions 
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of AO- and B-type operators. Remarkably, the ten different form factors involve only five 
independent linear combinations of the AO- and B-type coefficients as a consequence of 
helicity conservation of the strong interactions at leading power in the 1 / mf,-expansion. 
This implies 

-^V(E) = '^^^A,(E), UE) = ^UE) (23) 

mB + my 2E Ih 

up to power corrections [221 ? equality of the coefficients pertaining to pseu- 

doscalar mesons and longitudinally polarized vector mesons. The five pairs (Cj^°\ C^p^^) 
constitute the main result of the first matching step from QCD to SCETj. The 1-loop 
expressions can be inferred from ^Sj^]. They are collected in Appendix 1X1 



2.5 Physical form factor scheme 

Since the SCETj form factors E,a{,E) are not known, it has been suggested in [H] to define 
them in terms of three QCD (or "physical") form factors. Let 

fFF _ f cFF _ ^FF _ + my TTIb - my 

4p = /+, 4± = ^ , ^ 4|| = — — ^1 ^2, [^V 

TTIb + Tfly " LUi TUb 

which corresponds to [H] except for the longitudinal form factor The above definition, 
which has already been adopted in [25, is preferred, since it preserves the equality of 
the short- distance coefficients for the pseudoscalar meson and longitudinal vector meson 
form factors. 

In the physical form factor scheme we have again the two identities ((221); ^^e five 
remaining form factors read 

. _ p ^FF , (MBi) MBi) 

— /t = RriW + fC'ff^'' — C^^"^^ Rt) * Sp, 

rriB + rnp ^ •'^ ^ 



^n-T, = Rp + (Cjf ) - C|f ^' Rr) ^ S|, (25) 

Here the asterisk is a shorthand for the convolution integral over r. The ratios R of 
AO-coefficients take much simpler expressions than the individual coefficients given in 
Appendix |X1 Up to one loop [H] 



MAO) ^ 
^0 = 7^40) - ^ + 
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2-2£ 



Rt = 



C 



{AO) 
fr 



a 



(AO) 
Ti 



c 



(AO) 



1 + 



1 + 



ln4 + 2i 



In 



(26) 



with 



2E 



rrib 



2E mi, 



(27) 



Cp = 4/3, and u the renormahzation scale of the QCD tensor current. In the phys- 
ical form factor scheme there are only three non-trivial ratios R and three non-trivial 
combinations of B-coefficients. 



3 Jet-functions 

In this section we turn to the main part of this paper, the calculation of the SCETj form 
factors Sa(r, -E). Technically, this amounts to matching the B-type SCETi currents, 
<Jk{T)i whose matrix elements define the 'Ea{T,E), to four-fermion operators in SCETn. 
These four-fermion operators factorize into a product of two light-cone distribution am- 
plitudes resulting in the desired expression Q. That this can be done follows from [T2], 
where it has been shown that at leading power in the heavy quark expansion the B-type 
currents match only to four-fermion operators with convergent convolution integrals. In 
terms of operators, we derive the matching relation 

dudv Mt;v,u;) U^W,) (sn+) {W}C) (0)1 _ f (g.^) {tn^)^^5 {Y,^K) (0) 



2 «.\ C^/\ /J prp 



2 



FT 



+ . . . (28) 

to the 1-loop order, where the ellipses denote terms that have vanishing matrix elements 
between B mesons and pseudoscalar or vector mesons, or are power-suppressed in l/rnt, 
will be defined after ()57|1 . and the subscript "FT" denotes a Fourier transformation 
with respect to the light-cone variables s,t that will be made more precise later (see 
fl45|) ) . The functions Jk are the short-distance coefficients of the SCETn operators, which 
contain the hard-coUinear effects from the scale (m^A)^/^ integrated out in passing to 
SCETn. These short- distance coefficients will be called "jet-functions". 

3.1 Set-up of the calculation 

The jet-functions Jk{T\v,uj) are extracted by computing both sides of (j2HI) between 
appropriate quark states. We therefore consider the four-quark matrix element of the 
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operators Jk{T), 

A.(r;t;,^) = {q{p\)q{p'2)\Jk{T)\q{l)h{m,v)), (29) 

where k = 1,2,3, and the momenta /, p'^ {i = 1,2) are soft and colhnear, respectively. 
q denotes a hght quark, b the 6-quark. The quark and anti-quark in the final state 
may be of different flavour, but the flavours of the initial and flnal state anti-quark 
are identical. The quark-antiquark initial and flnal states are assumed to be in a colour- 
singlet state. Since the operators on the right-hand side of fj28p do not contain derivatives, 
the transverse momenta of the collinear states can be set to zero, and we can deflne 
p[ = vp' = vEn^, p'2 = vp' with v = 1 — v. Likewise for the soft states the momenta can 
be taken to be m^v for the heavy quark and / = uJUj^jl for the light anti-quark. The 
four functions Ha(r, E) deflned in (fTHj) . (plj) correspond to setting = (75)7^^ (Sp, Hy), 

Tfc = 757i^x7^^ (S±), and = 757''^7i'_l (-^)- 

The operators i7fc(r) generate momentum-space vertices with an arbitrary number of 
n^Ac gluons due to the Wilson lines Wc- Of these only the one- and two-gluon vertices 
are needed in the 1-loop calculation. The corresponding Feynman rules read (collinear 
quark lines dashed, all gluon momenta are out-going) 



k, fJ- Pi J n+k\ ( n'^k'l 




T^Tfc, (30) 



n^k 



k, p, A 



Ffc, p 




p'l 



ki, p, A k2, (T, B 



5{ T 



n+k2 
ri-^p' 



6 T 



?i+(fci + k2) 



n+k2 

n'+{ki+k2f^ 
n+{ki + k2) 



+ 5[t 



n+{ki + k2] 



{ki <^ k2, A <^ B, p <^ a). 



(31) 



In light-cone gauge n+Ac = the variable r corresponds to the fraction of total collinear 
longitudinal momentum n^p' carried by the transverse hard-coUinear gluon. 

In addition the calculation requires the collinear interactions from the leading-power 
SCET Lagrangian, 



C = U in_D + ilp^^ 



1 



-iiPi 



-e - - tr (FrF^.c) 



+ hJvDshy + qsilpsqs ~ 2 {FrFfius) 



(32) 
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as well as the sub-leading interaction [12] 



4^ = QsWliIpU - ^iV^c^cQs (33) 

that converts the soft spectator anti-quark in the B meson into an energetic, collinear 
anti-quark. The Feynman rules for the collinear interactions can be found in ^^1; while 
()33|) implies the vertices (collinear quark line dashed, soft quark line solid, gluon momenta 
outgoing) 

k, p, A 

^ — i,.T-(,l-^). (34) 

ki, p, A k2, a, B 




+ {ki ^ k2, A ^ B , p ^ a). (35) 

We note that n^p = n-^.{ki + /C2) and p± = {ki + /c2)±, since the corresponding soft 
momentum components are neglected in coUinear-soft interaction vertices (multipole 
expansion) . 

3.2 Unrenormalized matrix element 
3.2.1 Tree 

The tree contribution to (j29|) is shown in Figure d The gluon momentum (outgoing 
from the operator vertex) is given hy k = p'^ — I, and with k'^ = —vn^p'n_l = —2Eu!v, 
we obtain 

4°^ (r; ^) = Sir - t;) F, ® 7^ (36) 



mtv (s) ; p[ (c) 



(he) 



I is) p',ic) 



Figure 1: Tree diagram. 
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Figure 2: 1-loop diagrams. The first diagram summarizes all gluon propagator correc- 
tions. Diagrams without numbers are omitted when the calculation is done by expansion 
of QCD diagrams. 

where Fi ® F2 means Uc{p'i)TiUh{mbv) v{l)T2Vc{p'2)- The heavy quark spinor satisfies 
ifuhirribv) = Uh{mbv), and for the collinear spinors r^^vdp^) = Uc{Pi)'4- = 0- 

3.2.2 One loop 

A generic 1-loop diagram in SCETi contains contributions from the hard-coUinear, 
collinear and soft momentum region. For our external momentum configuration soft 
and collinear loop integrals are scaleless, so the diagram computation extracts the hard- 
coUinear contribution. The SCETj diagrams are shown in Figure |21 omitting diagrams 
that are obviously scaleless. We use dimensional regularization {d = A — 2e) for both 
ultraviolet and infrared singularities, hence scaleless integrals vanish. 

We computed the 1-loop diagrams in two different ways. First, we used the SCETi 
Feynman rules as given above and computed the diagrams as shown in the figure. Second, 
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we computed the matrix element (j^^ with full QCD Feynman rules, but expanded the 
Feynman integrand under the assumption that the loop momentum is hard-collinear, 
and the external momenta coUinear and soft, respectively. This method, also known as 
the "strategy of expanding by regions" [23 IIE] , gives precisely the same result as the 
first computation, but it turns out to be algebraically simpler, because it avoids having 
to use the more complicated vertex Feynman rules generated by the SCET Lagrangian. 
There are also fewer diagrams to compute. There are no two-gluon qqgg vertices in full 
QCD, and the unnumbered diagrams in Figure El are absent. In both computations we 
write ^ ^ ^ 

- — -7 = — -dn,kdn_k- — -7-77, (37) 

and first perform the n_/i;- integral by contour integration. The fc^-integral reduces to 
a conventional Feynman integral; the n+fc-integration is eliminated by the 5-function in 
(jnO)) or with the exception of diagrams such as (1), (5) or (6). 

It is convenient to perform the calculation without specifying the Dirac matrix of 
the SCETi operator. The unrenormalized matrix element (|29|) can be written as 

in terms of Dirac structures that cannot be reduced further in d dimensions. The notation 
is such that for any quantity the superscript (0) denotes the tree and (1) the 1-loop 
contribution. The coefficients of all three structures are infrared divergent, but only A^^^ 
and -B^^^ are ultraviolet divergent. It follows that all B-type currents can be renormalized 
with only two independent renormalization constants. Specifically, inserting the three 
Dirac structures (fT7|) . we obtain 

^S-) = [A+{d- 2)B] (75)7;., ®i^-+C (75)7''"7'"7^" ® 1,.1x^1p.. 
4"'^) = A (75)7.x7mx ® 7^^^ + (4 - d)B (75)7/.x7.^ ® 7'^^ 
+ C (75)7''^7^^7.x7''^ ® If^^lx^lp^, 

At^ = [A + {d- 2)B] (75)7mx7.x ^Y^+C {l5)Y^l'^Y^lu, ® 7mx7a,7px • (39) 

Here and in the following we use anti-commuting 75, and the bracket around 75 refers 
to the two cases, where 75 is or is not included in F^. From this it can be seen that the 
UV divergent parts have the same Dirac structure as the original operator. Hence, the 
operators J7fe(r) and all the B-type current operators do not mix under renormalization 
(in the basis adopted here). The renormalization constant for the operator J7fc(r) with 
Fi 3 is related to the divergent part of A^^^ + 2B^^\ the one for Jfc(r) with F2 to A^^\ 

3.3 Ultraviolet counterterms 

The counterterm diagrams are obtained from the tree diagram by insertions of a coun- 
terterm vertex into the gluon propagator, the quark gluon vertex and the operator vertex. 
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Finally, the on-shell matrix element must be multiplied by the propagator residue factors 
Rl/^Rq^^'^R^. They are related to the field renormalization constants by 

Rx = (40) 

where is the renormalization constant of field X in the on-shell scheme (by definition 
Rx^ = 1) and Zx is the field renormalization constant in the MS scheme. 

The Lagrangian counterterms are standard. As regards the operator counterterms, 
we note that the three operators j7/c(r) do not mix (see above), hence the renormalized 
operator is related to the bare operator (expressed in terms of bare fields) by 

JiAr) = jdr'Z\\{T,r')J^-{r'l 

Mr) = Jdr'Z^{r,r')Jt%T'), (41) 

which defines the operator renormalization kernels. Here we used that Ji{t) and Jz{t) 
renormalize identically. 

Putting together all the renormalization constants, the on-shell ultraviolet-renorma- 
lized matrix elements of Ji and J-^ follow from by the replacement 

A+{d-2)B ^ A+{d-2)B + Z\^^ ^ 

+ (2Z« + 1 ^) + + 2Zf (^))) AW, (42) 

while the renormalized matrix element of J2 is the second line of ()39|) with 

A + Z^^^A^^^ 

'2ZW + - ^) + ZOS(i) + 2Zf (^))) (43) 



The asterisk denotes convolution in the variable r' as in the definition of the operator 
renormalization kernels. The on-shell field renormalization factors are all equal to 1 
in dimensional regularization, so Zx^^^^ = 0. Zg is the standard MS strong coupling 
renormalization factor, g^^^"^ = ZgQs with {Ca = 3, Tj = 1/2) 

= 1 - /?o = ^ - -%T,. (44) 

^ Svre 3 3 ^ ^ ' 

The matrix elements (j39|l with the substitutions (j^2|) . ()43|1 are now ultraviolet finite, but 
infrared divergent. The infrared divergences are reproduced by the SCETn computation 
of the matrix element, resulting in a finite jet-function. This can be used to determine 
the operator renormalization kernels Z\\ and Z^ alternative to the direct computation 
of the operators' ultraviolet singularities performed in [201 • That is, after extracting the 
jet-function, we shall obtain the renormalization factors by requiring that they render 
the result finite. 
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3.4 Matching to SCETn and extraction of the jet function 
3.4.1 SCETn matrix element 

To obtain the jet-function, we need the SCETn matrix elements on the right hand 
side of the matching equation (^Hj) to the 1-loop order. In the absence of collinear- 
soft interactions in SCETn the four-quark operator factorizes into a coUinear and a soft 
bilinear. We define 



FT 



27r 



ds e~ 



Piuj) 



FT 



1 

and the quark matrix elements 



4~ 



(45) 



^bq{uj\uo) = {^\P{uj')\q{l)hM)- 



(46) 



The Dirac matrices will be determined later by the Fierz transformation of the first 
Dirac matrix product on the right-hand sides of (jSH)- Note that the hadronic ma- 
trix elements of Q and P are precisely the meson light-cone distribution amplitudes. 
Collinear-soft factorization'^ in SCETn means that the matrix element of Q\r^jJ\{v)P{uj), 
which appears in (PHj) . factorizes into 

{M{p')\Q[Tl]{v)P{u)\B,) = {M{p')\Q[Tl]iv)\0){0\P{uj)\B,), (47) 



which is a product of light-cone distribution amplitudes. 

With these definitions the tree quark matrix element of the SCETn four-quark oper- 
ator is given by 

(g(p'i)?(P2)lQK](^0^(^')k-(OKwO)^°^ = S{v - v')5{u - J) ® t^^, (48) 

where now ri®r2 means Uc{p'i)TiVc{p'2)v{l)V2Uh{jnbv). The "tensor products" ® and 
(8) are related by a Fierz transformation as discussed below. 

^We recall here that in general the apparent factorization of coUinear and soft degrees of freedom 
in SCETii is not valid ^1 due to the non-existence of a regulator that preserves factorization. 
The unregulated integrals correspond to endpoint-divergent convolution integrals. However, it has been 
shown that for the particular case of the matrix elements of the B-typc currents, the convolution integrals 
must be convergent ^5], so collinear-soft factorization is valid here. 
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The SCETn computation of the matrix elements is very simple, since in the 
coUinear sector SCETn is equivalent to full QCD, and in the soft sector to heavy quark 
effective theory The external momenta do not allow to form a non- vanishing kine- 
matic invariant, hence all loop integrals are scaleless and vanish. The matrix elements 
are given by tree diagrams including tree diagrams with counterterm insertions. We can 
therefore write 



uj) 



Zp{uj',uj) ^75, 



(49) 



with Zq{v',v) and Zp{llj',u!) the renormalization kernels that relate the renormalized 
operator to the bare operator expressed in terms of the bare fields, 

r-l 



dwZQ{v,w) Q[r^ 



c 1 bare / 



W 



P{uj)= du' Zp{uJ,Uj')P^'''''{ 

Jo 



UJ 



(50) 



The required 1-loop renormalization factors have been computed in other contexts. 
For the coUinear operator Zq{v,w) is the Brodsky-Lepage kernel 

Zq{v,w) = 6{v-w) + Z^q\v,w) + ..., 



zS\v,w) 



agCp 



2Txe 



1 

WW 



_ 9(w — v) _ 9(v — w) 
vw h vw 



w — V 



V — w 



- s{v - w) + A (—e{w -v) + —e{v - w^ 

2 \w w 



(51) 



where A = 1 applies to = (75)1 (pseudoscalar meson, longitudinally polarized vector 
meson distribution amplitudes) and A = to = 7^ (transversely polarized vector 
meson distribution amplitude), and the plus-distribution is defined for symmetric kernels 
/ as 



dw[f{v, w)]+ g{w) = J dwf{v, w) {g{w) - g{v)) . (52) 
Similarly, the renormalization of P{uj) has been worked out in [211 with the result 
Zp{uj,uj') = 6{uj - 00') + zP{uj,uj') + . . . , 

agCp 



Zp [UJ, UJ ) 



- + -\n^--)6{uj-uj') 
e UJ 2e 



2uj / 1 e{uj' -iu) 1 e{uj-uj'y 

e \u' uj' — UJ UJ UJ — uj' 



(53) 
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3.4.2 Extraction of the jet-function 

The jet-function is extracted from the quark matrix element of 



Ak{r;v,uj) 



duo' / dv' Jkir; to', v') <^bqi^', to) ^gqiv', v). 



This gives 



(54) 



(55) 



at tree level, and 



4'^ (r ; ^, ^) ^r,® ^75 = 4' V; ^, ^) 



dui' Zp^ {u, uj) jj^^ (r; u', v) 



(0), 



dv'Z^^\v',v)jf\r-u,v') 



(56) 



at the 1-loop order. Ji^\T\v,uj) is the uhraviolet-renormalized but infrared divergent 
1-loop matrix element of JkiT). The subtraction on the right-hand side precisely cancels 
the infrared divergences, so that the short-distance jet-function is finite as it should 
be. There is however a difficulty in executing the subtraction in the form of ()56|) . since 
J^^\t]v,uj) is expressed in terms of the spinor product ®, corresponding to SCETn 
operators with spinor indices contracted in [^h^] [qsS]-, while the left-hand side and the 
subtraction term involve the product ®, corresponding to operators with spinor indices 
contracted in [^^] [qgh^]. The standard Fierz identities that relate these structures are 
valid only in four dimensions. We shall discuss this issue together with the reduction of 
evanescent Dirac structures appearing in in the following subsection. 

None of this applies to the tree-level equation where all quantities are finite. 
We can therefore apply the four-dimensional Fierz identities to the tree-level terms in 
((23), so that for 

k = l: (75)7/.x ® 1^"- = (±1) y75(75) ® + • • • , 

= 2 : (75)7.x7/.x ® 7''^ = -y7.x 75(75) ® ^75 + • • • , 



3 : 



(75)7Mx7i^i ® 1^ 



(57) 



The Fierz transformation produces a number of 4-fermion operators with different Dirac 
structures. In the following we discuss only those operators (jet-functions) that con- 
tribute to the decay of a pseudoscalar B meson. The ellipses denote terms involving /i-/2 
and 2 71/^^(75) to the right of (8>, which vanish when the matrix element (0|g<j[- • ]hv\B) 
with the pseudoscalar B meson state is evaluated. Hence, these terms will not be consid- 
ered further. The upper (lower) sign refers to the operators on the left-hand side without 
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(with) the 75 factor. Here and below the Fierz identities are given for four-fermion op- 
erators and the extra minus sign from the field permutation relative to the identities for 
matrix products is already included. Comparison of (jFfj) with the definition (jl^ of Q[r^] 
determines the Dirac matrix in the collinear SCETn operator. For k = 1, we have 
Tfc = (75)7mx and the corresponding = (±1)75(75); for k = 2,Tk = (75)7i^x7mx and 
F^ = —71^^75(75); for fc = 3, Ffc = (75)7^^7,^^ and F^ = 0, i.e. there is no contribution 
for pseudoscalar B mesons. Hence, comparing (|Ffj) to (IHBj) leads to the tree-level 
jet function 

4»'(r;.„..,^-f£5^*(.-.) (58) 

for k = 1,2, and for k = 3. 

Using the 1-loop expressions for the SCETn renormalization kernels, we can evaluate 
the subtraction term in square brackets in ()56|) with the result 



Nc 2Eujv 



3.4.3 Evanescent operators and Fierz transformation 

We now discuss the reduction of the Dirac structure and the Fierz transformation. We 
consider in detail the case F^ = 7^^. According to the first equation of (jH^ the UV 
renormalized matrix element of the corresponding SCETi current J7i(r) is 

Aiir-v,uj) = A-f^^^Y^ +CY^1^^Y^ ^If^^lx^lp^, (60) 

where A stands for the right hand side of (j42p . The second Dirac structure reduces to a 
multiple of the first one in four dimensions, and the first one satisfies the Fierz relation 
fl57|) in four dimensions. In the following we discuss separately the reduction of the 
Dirac structure and the Fierz transformation in d dimensions, although it is possible to 
perform both reductions in a single step. 

Indicating only the field content and Dirac structure we define the operators Oi = 

{^lti±h^){QsY^O and 0[ = {^Y^'y^^Y^K){qslf,^lx^lp^O^ and rewrite the previous 
equation as 

Aiir;v,uj)= (Oi)(°) + CW(0;)(°), (61) 

where the tree-level matrix elements of O^'^ just reproduce the Dirac structures, i.e. 
(Oi)^^-' = UciPi)Y^Uh{mbv) v{l)'yfj_^Vc{p'2) etc. In d = A, 0[ = 40i, so we define the 
evanescent operator E = 0[ — f{e)Oi, where /(O) = 4, but otherwise /(e) is arbitrary. 
Hence, 

Ai{t-v,u) = + + /(e)C«] (Oi)(°) (62) 

Since C^^^ has a 1/e infrared divergence, the coefficient of the "physical" operator Oi 
depends on the up to now arbitrary prescription /(e) in the definition of the evanescent 



20 



operator. The jet-functions are obtained by expressing this equation in terms of the 
renormahzed operator matrix elements. To the 1-loop order it is sufficient to 



use 



(63) 



which results in 



Ai(t; v,uj) 



i(0) + id) + /(e)C(i) - i^Mgl (Oi) 



{E) 



J, 



(0) 



Oi 



(1) 



Oi 



(Oi). 



(64) 



The coefficients of the operator matrix elements in the first line of this equation are now 
infrared-finite short-distance quantities, and the equation is interpreted as an operator 
matching relation valid also when the matrix elements are taken between hadronic final 
states. Hence we can take the limit d — > 4 in which (E) = 0, since the tree-level 
matrix element vanishes in four dimensions.^ The 1-loop jet-function Jq^ can be read 
off from It depends on the choice of the evanescent operator, i.e. the order-e 

term of /(e), but so does (Oi) as can be seen from (jU^ . and the physical amplitude is 
scheme-independent. 

The jet-function Jo^ is not the desired result, because instead of Oi we must use the 
Fierz-transformed operator 



Pi 



Qs—l^K 



(65) 



where the ellipses denote terms that are irrelevant for us, as in fj57|) . This Fierz-ordering 
is uniquely singled out by collinear-soft factorization in SCETn- Radiative corrections 
to Pi occur only within the collinear factor or within the soft factor, hence the Dirac 
structures in arbitrary loop diagrams can be reduced to the tree structure Pi. It follows 
that whatever evanescent operator one may write down in this Fierz-ordering decouples 
from Pi, and can simply be ignored. Hence, even in d dimensions we have 



t(0) 
•J p. 



(1) 



{Pi). 



(66) 



with Jp^ the desired 1-loop jet-function. Comparing this equation to we have 
= Jq^i if the renormahzed matrix elements (Pi) and (Oi) are equal. In general this 
is not the case, because the Fierz identity that relates Oi and Pi in four dimensions is 
not valid for d ^ A. Now, the renormalization scheme for Pi is fixed by the requirement 
that the collinear and soft matrix elements coincide with the standard MS definition 
of the light meson and heavy meson light-cone distribution amplitudes, respectively, so 
(Pi) is completely defined. However, by choosing /(e), we can adjust the definition of 
Oi, such that the infrared-finite matrix elements are equal in the limit d — > 4. 



"^In higher orders one usually uses a non-minimal subtraction scheme to ensure that {E) 
this is not relevant to the present calculation. 



0, but 
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Since the jet-functions are independent of the infrared regularization, any one can 
be chosen for the calculation, and it is convenient to compare the matrix elements of 
Oi and Pi by assuming an off-shell infrared regularization. The ultraviolet-renormalized 
1-loop matrix elements are 

(Oi) = (l + MS'f) (Oi)W+MS^^(i?)W, 

(Pi) = (l + M°f)(Pi)W. (67) 

There are no 1/e poles in these equations due to the use of off-shell IR regularization, 
so the evanescent term drops out for d ^ A. Furthermore, {Oi)^^'^ = (Pi)*-"^ by the four- 
dimensional Fierz-equivalence of Oi and Pi, hence we need to define Oi such that the 
difference Mq^ — Mp^ = 0. At the 1-loop order one finds that the self-energy corrections 
and the vertex correction, where the gluon is exchanged between the soft light quark and 
the heavy quark, drop out, because the coupling to the heavy quark is proportional to 
v'^ and does not introduce a new Dirac structure. Only the gluon exchange between the 
coUinear ^-fields can give a non-zero difference. In the case of Oi the part of this diagram 
that does not cancel in the difference Mq^ — involves '-fP^'y^^'yt^^ ® 7/ij^7Ax7p^, and 
its contribution to Mq^ is therefore proportional to /(e). In the case of Pi, the Dirac 
structure is 

7"-7^-^7/.,7a,®^ = id- 2)2(Pi)W. (68) 

Including the overall coefficient, we find 

M^^l - M°f = (/(e) -{d- 2y) (Pf) (69) 

where [C^^-']div is the divergent part of C^^\ In order for this to vanish, we must choose 
/(e) = (d — 2)^. In other words, we must define 

Y^^^^Y"- ® 7mx7a,7px = id- 2)2 7^^^ ® 7^^. (70) 
With this Jp^ = Joi, so returning to we have 

jff = \im (i(^) + {d- 2)2^(1^ - i^M},;)) . (71) 

The term A^^^Mp^ is nothing but the subtraction term ()59p that appears in ()56p. while 
A^^^ and C^^^ can be read off from the 1-loop calculation that leads to and (HH), so 
the previous equation gives the final result for the correctly renormalized and subtracted 
jet-function. The above argument can be repeated for = 757^^^, and one finds that 
the jet-function is identical to the one for F^ = 7^^ as was expected. 

The case of B-type currents 1/2,3 (r) with one uncontracted transverse index is slightly 
more complicated than the scalar case, because there are two physical and two evanes- 
cent operators. Corresponding to = {''i^)^^!^^ and (75)7^^-^7^-^ we define the two 
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physical operators O2 = {^{l5hu^lfiMilsl^^0 and O3 = {^{^5h^l^luMi1sl''^0■ 
Requiring that the renormahzed matrix elements of 02,3 equal the renormalized matrix 
elements of the corresponding operators P2,3 in the other Fierz-ordering fixes uniquely 
the prescription for reducing the Dirac structures multiplying C in (j39j) and ensures that 
the SCETn operators are correctly minimally subtracted. A short calculation analogous 
to the one discussed above gives 

(75)7''^7^^7.^7''^ ® 7mx7a^7px = (d ' 4)' (75)7.x7/.x ® 7^"^, 

(75)7'^7^^7''^7.x ® 7mx7Ax7px = (d - 2)' (75)7/.x7.x ® l"^- (72) 

Referring to (j37j) we see that the matrix element of O3 vanishes for a pseudoscalar B 
meson. Hence the B-type operator J^s^t) with Tk = {'l5)'lfi±1uj_ has a vanishing jet- 
function just as at tree level. On the other hand, inserting the first equation of (f?^ into 
the second of we find that the B-type operator (r) matches to JP2P2 with 

jg) = lim + {d- 4)2C(^) - A^m},;)) . (73) 

Here A^^^ and C^^^ are defined by (jH^ and PHj) and the subtraction term is given by 
fl56|) except that now A = must be used in the Brodsky-Lepage kernel. Since C-^-* has 
only a single 1/e pole, the term [d — 4)^C^^'' does not contribute to the final result for 
the jet-function. 

3.5 Final results for the jet-functions and B-type current re- 
normalization kernels 

We shall now denote the jet-function Jp^ {JP2) that arises in the matching of J7i(r) 
{J2{'t)) as J|| {Ja_). To the 1-loop order we write 

with a = II, _L. Applying the subtraction procedure described in the previous subsections 
we obtain the ultraviolet and infrared finite 1-loop corrections ja(r; v,a;). The resulting 
expressions still depend on the yet undetermined renormalization factors Z'^^ of the B- 
type currents through (jl^ . ()43|) . They can now be determined by the condition that 
jair-^v^uj) must not contain 1/e poles. We choose the MS scheme to be consistent with 
the definition of the B-type short- distance coefficients in |18j . 

3.5.1 Renormalization kernels 

We expand the Z-factors in pT|l as 

Z,(r,rO = 5(r-r') + ^^«(r,r'), (75) 
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and obtain 



with 




The anomalous dimensions of the B-type operators are derived from the renormahzation 
kernels in the standard way (see (jHH) below). Our result is in complete agreement 
with [20], where the anomalous dimension has been obtained by extracting directly the 
ultraviolet divergent parts of the SCETj diagrams.^ 



3.5.2 Jet-functions 

The 1-loop jet-functions read 

Jll (r; v,u) = A 5{r -v)+{Cf-^] [2Bl 



div -t)—[L + \ti 
vv * 



[V — TjT V[V — T] 



TT 



2vf 



L + In rr 



6{t-v) 



2{v - rf 



VVT 



L + ln 



+ 9{v-t) 



2{v - t) 



vv 



L + ln(f - r) 



(T — V)T 



TT 



VT 



IV - r 



In^- 



{v — t){v — t) t f — r 



vt\' 
vtJ. 



T — v) —IL + ln(r — v) + 



, V T 

— m — 

T — V T V 



^The variable u in |2U| corresponds to our t = 1 — r, their v to our r' . 



(78) 
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jx(r; v,uj) = A5{t-v)+[Cf-^] [2B]_ 



v\ V 



2f 



[T - V 



Op — 



eir -v)—[L + \n 

VT ^ 



VT 

T — v)f V — T 



V 

V 1 



VT 



2 / T 

+ 9(v-t)-(L + \n(v - r) + In 

V \ V — T r V 



+ &{T-v)—[L + ln(T-v) 

VT 



VT ^ V 

m 1 



r — V r 
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with 



A = Cf 



B 



[L + lxivj^ -"^[L + lnv 
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L + ln(r — v] 



9{v - t) 



L + \ia{v — t] 



2n^ 152 
~3 9~ 



(80) 



T — V \ J V — T 

and L = ln(n+p'n_///i^) = ln(2_Eu;//i^). The SU(3) group factors are Cp = 4/3, Ca = 3, 
Tf = 1/2, and nj denotes the number of hght quark flavours. Once again we find 
agreement with [T^ I20j . 



3.5.3 Hard-scattering form factors 

Here we express the hard-scattering (SCETi) form factors E{r,E) defined in (fTH|) . (pUjl 
as convolutions of the above jet-functions and light-cone distribution amplitudes. 

The light-cone distribution amplitudes of the light mesons follow from the matrix 
element of (5[r^](w) defined in (j^ . 



{M{p')\Q[n]{v)\o) 



-ifpE(t>p{v) 

. mye* -v 
-ifv\\E (t)v\\[v) 

-ifv±E (e*" - e* ■vn'^) ^v±{v) 



■75 



(81) 



The three cases correspond to M being a pseudoscalar meson, longitudinally polarized 
or transversely polarized vector meson, respectively. Similarly, the B meson distribution 
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amplitude is related to the matrix element of P{uj) in (jlSj) such that 

(O|PH|5„) = ^^0B+M (82) 

with Jb the HQET B meson decay constant (but defined such that is has mass dimension 
l)^ that is 

fB = K{^.)fB{f^) = ^^^^^^^^ (83) 
with -Fstat(At) the m^-independent decay constant in the static limit (HQET) and 

if(ri^l + ^(3.„!^-2) (84) 

a short-distance coefficient [30] . 
We can now evaluate 

(P(p')|Ji(r)|B.) = 2eI ^e-''^^^ {Pip'mwmiW^^^^,W,){rn+)KmB.) 
= J dudv J\\{T;v,uj){P{p')\Q[^-f5\ {v)P{uj)\B,) 

^ dujdv J\\{t]v,uj) fB<pB+{'^) fp<pp{v), (85) 

where we have used fl2Hj) and (jlTj) . Comparing this to the definition (fTHj) and (fTUj) gives 
the expression for Sp. Proceeding in the same way for the other form factors we obtain 

2p(r, E) = ^-^ * /p0p * 

— -II (^> ^) = T-^ /p0B+ * /y||0v|| * ^11 , 
my 4mb 

2±(r, E) = ^-^ * fvL<Pv± * 

Hx(r,E) = 0. (86) 

The asterisk stands for the convolutions in u and v as in (j85|) . This together with 
the explicit expressions for the 1-loop jet-functions is the main technical result of this 
paper. Using this result in (^3)) allows us to investigate numerically the corrections to the 
symmetry relations for heavy-to-light form factors at the 1-loop order. The subsequent 
sections are devoted to this numerical investigation. 

^With the conventions of heavy-quark effective theory our \By) corresponds to ^tob|B„), and the 
right-hand side of H82|l is ms-independcnt when expressed in terms of i^stat ■ 
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4 Numerical analysis 
4.1 Jet-functions 

It will be seen below that the jet-function appears in the form factors in the form of the 
integral 




which is normalized to 1 in the absence of the a^-correction. We now evaluate the 1-loop 
correction. 



4.1.1 Light-cone distribution amplitudes 



The light-cone distribution amplitude (LCDA) of the light meson, (puiy), is convention- 
ally expanded into the eigenf unctions of the 1-loop renormalization kernel, 



(88) 



where a*^ and C^^'^\v) are the Gegenbauer moments and polynomials, respectively. We 
define the quantity 



M 



1 dv 



V 



3 1 + E«n, 



M 



(89) 



n=l 



In practice the expansion will be truncated after the second term, since it is believed 
that the higher Gegenbauer moments are negligible, or accounted for approximately by 
phenomenological determinations of the first two moments. The LCDAs and Gegen- 
bauer moments are scale- and scheme- dependent. Our computation of the jet-function 
corresponds to the modified minimal subtraction (MS) scheme and the renormalization 
scale /i of the LCDA (not indicated by its arguments) is equal to the scale n that appears 
in the expressions (f7S|) . (fTHj) for the jet-functions. In particular the scale-dependence of 
{v-'^)m is given by 



an 



\nv r 



M 




vA 



(90) 



which follows from ()5H).'^ We recall that A = for transversely polarized vector mesons 
M, and A = 1 for pseudoscalar mesons or longitudinally polarized vector mesons. 

'''Note that (|5T|l describes the scale dependence of /a/(/'a/(w), and 

asCr A - 1 



d 



M- 
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The first inverse moment of tlie LCDA of tlie B meson, 



1 f'^ did , , , . 

— ^ / — (91) 

is a key quantity in exclusive B decays We define the averages 

{f{u))^\B —ct^B+{oj)f{u). (92) 

JO OJ 

The LCDA and Ab are scale-dependent. Our computation of the jet-function corresponds 
to the modified minimal subtraction scheme and the renormalization scale /i is equal to 
the scale /i that appears in the expressions for the jet-functions. The scale-dependence 
of 1/ Ab is given by 

which follows from (j53|l . The first term in the bracket comes from the scale-dependence 
of the static decay constant /b. 

Since uj is of order A, only logarithmic modifications of the first inverse moment 
appear at leading order in the l/mf,-expansion. It can be seen from (fTSj) and (f?^ that the 
1-loop calculation involves the two logarithmic moments (L), (L^) with L = \n{2Euj / n'^). 
The entire energy and scale-dependence of the 1-loop jet-functions is contained in these 
two quantities. We adopt a simple one-parameter model for the shape of the distribution 
amplitude p?T] . 

0B+(c^) = ^e-A-, (94) 

^B 

which relates the two logarithmic moments to the parameter Ab, 




L =ln ^ , {L') = W ^ + - 95 



with = 0.577216 . . .. 

The functional form ()94|) and the Gegenbauer moments are assumed to be given 
at some reference scale /io of order (m^A)^/^. This avoids having to evolve the me- 
son parameters from the hadronic scale A to the hard-collinear scale with the SCETn 
renormalization group equations.* 



4.1.2 Integrated jet-functions 

We proceed to the evaluation of la- The integrals dr ja{T]v,uj) are given in analytic 
form in Appendix El Integration over uj introduces the logarithmic moments (L), (L^) 
defined above. The final ^-integration can be done numerically, or term by term in the 

®The corresponding anomalous dimensions are given by H51() . H53() . 
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Gegenbauer expansion of the light meson hght-cone distribution amphtude. Up to the 
second moment, we obtain 



1 + 



M 



1 + af + {L') - 5.24 + 8.93af + 10.86af (L) 



M 



+ 



3.99 + 8.67af + 13.47af 



asjfJ') 3_ 

Air 



M 



{L') - 4.90 + 8.93af + 10.814" (L) 



+ 



0.73 + 6.78af + 11.48af 



(96) 



with^ Uf = 4, Tf = 1/2, Cp = 4/3 and Ca = 3. The analytic expressions of these 
convolutions are also given in Appendix FBI 

These results do not contain large logarithms, when /i is of order of the hard-collinear 
scale (mbA)^/^ ^ 1.5 GeV, i.e. (L), (L^) are of order 1 for E ~ rrih and fx ~ (mfoA)^/^. 
Since as(1.5 GeV)/(47r) is approximately 0.029, (L^) 2.5 and (L) ^ —1 (for typical 
parameters), the perturbative corrections to the jet-functions are about (20-50)%, de- 
pending on a = II , -L and the precise values of the Gegenbauer moments. We may therefore 
conclude that perturbative corrections to hard spectator-scattering are non- negligible. 
At the same time there is no sign that the series expansion is not well-behaved despite the 
comparatively low scale, lending support to the possibility of performing perturbative 
factorization at the hard-collinear scale. This is an important result, already mentioned 
in [201 ESj, since theoretical calculations of exclusive B decays in general rely on this 
possibility. The present calculation and the one in [TTH OH] are the first computations of 
quantum corrections to spectator-scattering. 



4.2 Renormalization group improvement of the C^^^^ coeffi- 
cients 

The complete hard-scattering term involves Jq dr [C^^ * J] , so we now turn to the evalua- 
tion of the C^^. With n of order (m;,A)^/^, the jet-function is free from large logarithms, 
but C^^ involves up to two powers of logarithms ln(mf,/A) per loop from the ratio of the 
hard to the hard-collinear scale. In the following we derive an expression that sums the 
leading-logarithmic and double-logarithmic terms to all orders in perturbation theory. 
We shall refer to this simply as the leading-logarithmic approximation (LL).^° 

^We assume four massless quark flavours throughout this numerical analysis for simplicity. This is 
not a good approximation for the charm quark, whose mass is of the same order as the hard-collinear 
scale. A more precise treatment would keep the charm quark mass in the fermion loop correction to the 
jet-function. This could be done without difficulty if such precision were required. 

"'^''In the literature on Sudakov resummation the analogous approximation is usually called "next- 
to-leading-logarithmic approximation" . We prefer the term "leading-logarithmic" , since, as in other 



29 



4.2.1 Solution to the renormalization group equation 

The renormahzation-group formahsm that accomphshes this summation is standard and 
has been applied to the B-type SCETj currents in [SO]- We shall recapitulate the relevant 
expressions to define the notation. In the following we drop the superscript "Bl" on C^^^^ 
and denote by C{E,T;fi) a generic short-distance coefficient. With the renormalization 
kernels ()76p. (f77j) the renormalization group equation is derived from the requirement 
that C (E , t; n) J'^^^^jj) is independent of the QCD/SCETi factorization scale. This 
implies 

r; fi) = -r,,3p(«s) In ^ C{E, r; /i) + /' rfr'7a(r , r) C{E, t'; fi) (97) 
dfi 2E Jo 

with 

rcuspK) = Er„(^j (98) 

the so-called universal cusp anomalous dimension. Comparison with ()76p . ()77p gives 
To = 4:Cf and 

7a(r,r') = -^[z,{T,r') + A,z,{T,T')], (99) 

where = 1 for a = \\ and Aq = for a = _L. Due to the presence of double logarithms 
we also need the two-loop cusp anomalous dimension 



«^ ^'1C.-?„,T,V ,100) 



9 
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The short-distance coefficients of the B-type operators (jH)) to |TT]|) and the coefficients 
appearing in (j2^, ()22|1 evolve with the anomalous dimensions 7a(T',r) as follows: 



^{Bl) ^{51)1-3 ^{Bl)l,2,5-7 ^(Bl) ^(Bl) ^(Bl) 



a = ± 45^', C?'^'-', C{f'\ 4f ^- (101) 
The general solution to the renormalization group equation (|97|) reads 

C{E, r; /i) = e-^^^^^-'^) dr' Ua{r, r'; /x) C(E, r'; /x,), (102) 



where 



S{E;fih,M = / das—Tz — ^ / 

JaJuh) P(a_.) J a 



j„ rcusp(a«) r= c?a^ , , /ih, , rcusp(«. 



Tr7 / + / TV (103) 

P{as) Jasit^h) PW 2E Jasif^h) p[as) 

applications of renormalization-group improved perturbation theory, the approximation requires only 
the 1-loop anomalous dimension of a generic operator. The complication from double-logarithms is 
reflected by the presence of the so-called cusp anomalous dimension, for which the 2-loop coefficient is 
needed already in the leading-logarithmic approximation. 
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with the QCD /5-f unction 



n+l 



/9o 



11 



-Ca - i^rifTf, pi 
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The evolution kernel Ua{T, r'; fih, f^) satisfies the integro-differential equation 

/i3-f^a(i", t'; fih, /i) = / dT"-ia{r", t) Uair", t'; fih, fi) , 
da Jo 



(104) 



(105) 



with initial condition Ua{T, r'; fit, jJ^h) = S{t — t'). To sum the large logarithms the initial 
scale Hh should be of order rrib, and the evolution ends at fi of order (m;,A)^/^. 

Several simplifications occur in the leading-logarithmic approximation. Eq. (|103p can 
be integrated to 
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(106) 



with r = as(/i)/as(/i/i) > 1. Furthermore, the initial condition is given by the tree-level 
expression for C{E,t; ^h), which is independent of r (and fih), hence the r' integration 
in ()102j) can be done. We then have 

C(^^)(E,r;/i) = e-^(^^'^-'^)f/,(r;/i,,/i)CW(i?), (107) 

where C'^'^^E) is the tree coefficient, and Ua{T; fih, fJ-) = lo dr' Ua{T,T'; fJ') satisfies 

d 

A*3-f4(T;/ife,/i) = / dr' -fa{r',T) Uair'] fih, fJ') , (108) 
d/i Jo 

with initial condition Uaij] Hh, Hh) = 1- This equation must be solved numerically. We 
always use two-loop running of ois{ij), and put /i/j = = 4.8 GeV. As input we take 
as(4.8GeV) = 0.215, which gives 0^(1.5 GeV) = 0.359 (four massless fiavours). The 
result of this integration is shown in Figure El for /i = 1.5 GeV. We have found that the 
solution to 



given by 



/.-^t/:PP(r,/.,,/i) 
dji 



[/r(r,/i,,/x) 



dT'-fa{r',T) 



U:^^{T,fih,fi), 
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-7a(r)/(2/3o) 



(109) 



(110) 




Figure 3: Ua{T, fih, fJ') ior fih = nib = 4:.8 GeV and fj, = 1.5 GeV. The upper curves refer 
to U\\, the lower ones to U±. Sohd hues: exact numerical integration. Dashed lines: 
approximate solutions. 

with ^7a(r) = /o dr'7a(r',r) and 

7||(r) = -Cf + 4:{Cf- ^] , 



2 / r 

7±(r) = -C^f^— + lj +4(^C^-^j Inf + ^j (111) 

provides a very good approximation (better than 1%) to the exact solution, provided 
one uses 1-loop running for with as(4.8GeV) = 0.215 in the approximate solution. 
The approximate expressions are also shown in Figure |S1 



4.2.2 NLO+LL approximation 

We are now in the position to give expressions for the B-type short-distance coeffi- 
cients C{E,t; fi), which include the complete 1-loop correction as well as the leading- 
logarithmic terms. The formula is 



+ C^'\E) [e-'Ua-l] (i?,r;/i). 



{E,T;fx) 



(112) 



The meaning of the four terms on the right-hand side is as follows: the first and second 
terms are the tree and 1-loop coefficients, respectively. Together they constitute the 
next-to-leading order (NLO) approximation to C{E,t; fi). The fourth term is the sum 
of leading-logarithmic terms to all orders minus the tree. Finally, the third term subtracts 
the logarithmic terms already included in the full 1-loop correction C^^\E,t; fx). The 
subtraction is given by 



(E,r;/i) 



Atx 



2 



In^ f^Usinf^Vn {^^ 
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+ .n(A)-,„Mj (113) 

To analyze the structure of the correction, we display in Figure |3] the following ap- 
proximations, all normalized to the tree coefficient: (i) tree plus the logarithmic terms 
at 1-loop (dash-dotted); (ii) previous approximation plus the non- logarithmic term, i.e. 
the complete next-to- leading order result (dashed); (iii) previous approximation plus 
the sum of leading- logarithmic terms at order and beyond (solid). In the numerical 
implementation we set Hh = rrih = 4.8 GeV, E = mf,x/2, and regard the coefficients as 
functions of energy fraction x and the convolution variable r. Since E must be of order 
rrib, X cannot be chosen too small. We take x = 1 and x = 0.6 as representative examples 
of large (maximal) and small energy of the light meson, and evolve to /i = 1.5 GeV. We 
also fix the scale of the QCD tensor current operator to u = rriiy = 4.8 GeV. Figure |3 
shows four of the five combinations, Cx , ^ G {/+, /o, /t, Ti, V}, which appear in the 
hard-scattering contribution to the vector and tensor current form factors (PT|) . ((221) • 
Only C\P^^ is not shown, because its tree coefficient vanishes, hence only C^^\E, r; /i) in 
()112j) is non-zero. 

The following observations can be made from the figure: a) the logarithmic term at 
order (dash-dotted lines) is a very poor approximation to the full coefficient (dashed 
lines), especially for c!j^^\ which is the only one of the four coefficients shown involving 
the transverse anomalous dimension 7_|_. b) Except near the endpoints r = 0, 1, where 
the relative correction diverges, the typical next-to-leading order correction from the 
hard scale mt is of order 30%. The endpoint singularities are logarithmic and disappear 
when the correction is folded with the light-meson distribution amplitude (integration 
over r). c) The effect of the logarithmically enhanced terms beyond the order as is 
neghgible (difference between the solid and dashed lines). It is largest for Cj^^ towards 
larger r since here U±{t; fih, f^) is significantly different from 1, see Figure El 

4.3 Spectator-scattering correction 

According to (EH), ((221), and ((SH|) the form factors Fx{E) are given by 

Fx{E)=&f^UE) + Hx{E) (114) 
with the spectator-scattering term 

Hx{E) = /b/m / dv(j)M{v) / duj(j)B+{Lj) 

Arrih Jo Jo 

X f'dTd^'\E,T)Ja{T-V,u). (115) 

Jo 

Here a = || (or P in case of C,a) for X = {/+, /q, /t} and a = ± for X = {V, Ti}. We have 
now assembled all the pieces required for the evaluation of Hx{E) at order (1-loop). 
From now on we set mB/mb in ()115|1 to 1, since the difference between the meson and 
the quark mass is a power correction beyond the accuracy of the present calculation. 
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Figure 4: The short- distance coefficients Cx {E,T,fi) relevant to the form factors X 
at two representative energy values (x = 1 (left) and x = 0.6 (right)) at /i = 1.5 GeV 
normalized to the tree approximation. Dash-dotted: logarithmic terms at order 
dashed: full NLO approximation; solid: NLO plus logarithmic summation. 
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At the leading order we insert the tree expressions for the B-type coefficient function 
(hereafter we again drop the superscript "Bl") and the jet-function and obtain 



HPiE) 



Cf{E). 



2E\ 



B 



(116) 



Here as before the superscript "(0)" refers to the tree approximation. This agrees with 
the results of jHj- 

To obtain the next-to-leading order result including the renormalization-group sum- 
mation, we insert ()112j) and the jet-function into ()115j) . and neglect cross terms of order 
in the product Cx * Ja- The result is 



Hx{E) = Hf{E) X 
1 
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{E,v) 
{E, v) 



;ii7) 



The second term in the bracket is the 1-loop hard correction; the third comes from the 
jet-function and is defined in (j^Bj) : the fourth and fifth are related to the renormaliza- 
tion group summation as in pi2j) . The integration of the second term can be done 
analytically, but the expressions are lengthy. They are given for selected short-distance 
coefficients and for the integration with the asymptotic distribution amplitude in Ap- 
pendix El It is as straightforward to perform the integration over 0m (w) numerically. 
The integration of the subtraction term is elementary and given by ()113|) together with 
(jlllj) . The integration of the last term can only be done numerically using the numerical 
solution of the integro-differential equation for Since S is independent of w, one needs 
the integrals 



^ dv 

— 0A/(y,/i) U\\{V, flh, fJ-) 
V 

— (puiv, n) U±{v, Hh, IJ) 

V 



3.037 + 3.058 af + 3.051 af + 



2.795 + 2.980 af + 3.003 + ... . 



:il8) 



The numerical values are given for /i^ = 4.8 GeV and /i = 1.5 GeV. 

To illustrate these results we consider the three coefficients relevant to the form 
factors in the physical form factor scheme defined in (j25|) . Let 



(Bl) 



0+ {r,E) 
^ r,E) 
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(Bl) 



4fy^(r,i?) 



(j{Bl)f^ _ ^(Bl) 



Ci'''\T,E)Ro{E), 
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(119) 
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Figure 5: Hx{E)/Hx \E) for X = 0+,T+,TiV. The two upper curves represent the 
complete next-to-leading order result with (solid) and without (dashed) the renormal- 
ization group summation. The middle (dash-dotted) line shows the jet-function (hard- 
coUinear) correction alone, the lower (dashed) line the hard correction alone. 

Choosing = 0.35 GeV, x = 0.85 (corresponding to light meson energy E = xmB/2 = 
2.24 GeV or momentum transfer = 4.18 GeV), and asymptotic distribution ampli- 
tudes, the curly bracket in ()117|) evaluates to 

X = 0+ : 1 + 0.283 [C] + 0.371 [jet] + 0.059 [log, a,] - 0.073 [all logs], 
X = T+: 1 + 0.213 [C] + 0.371 [jet] + 0.059 [log, a,] - 0.073 [all logs], 

X = TiV : 1 + 0.209 [C] + 0.268 [jet] + 0.169 [log, a,] - 0.147 [all logs], (120) 

where the five terms correspond to the five terms on the right-hand side of ()117p . We 
observe that the hard correction [C] and the jet-function correction [jet] are of similar 
size, while the sum of higher-order logarithms (the sum of the last two terms) is at least a 
factor of 10 smaller. The total correction to the tree result amounts to an enhancement 
of (50-70)% of the spectator-scattering effect. These features are independent of the 
value of E as can be seen from Figure which displays the weak energy-dependence of 
the spectator-scattering correction normalized to the tree result. 

The dependence of these results on the hadronic input parameters A^, af^, is 
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roughly as follows. As enters the relative correction through the moments (j^J^ and there- 
fore affects the jet-function terms only. Choosing = 0.5 GeV (0.25 GeV) changes the 
number 0.371 to 0.295 (0.452), and 0.268 to 0.195 (0.346), an uncertainty characteristic 
for all energy fractions x. Furthermore, there is an uncertainty due to the model for the 
shape of the distribution amplitude that correlates the logarithmic moments with Ag, 
which we do not attempt to quantify. There is a larger dependence of the tree result 
Hx\e) on Xb, since it is inversely proportional to A^. Positive Gegenbauer moments 
increase the tree result and the relative next-to-leading order correction. This can be 
seen from (j96|) for the jet-function correction. For = 0.2, the total relative next-to- 
leading order correction increases from 64% for X = +0 (50% for X = TiV) to 75% 
(62%). Finally, there is a dependence on the renormalization scale /i, which we fixed to 
1.5 GeV. In order to estimate this dependence, one must fix the hadronic input parame- 
ters at some /iq and evolve them to /i using (j^ . Since the scale-dependence of the 
hadronic parameters is within their uncertainty, we do not perform this estimate here. 

5 B decay phenomenology 

In this section we discuss three applications of our results to B decays. We restrict 
ourselves to decays to pions or p mesons, since the results for kaons are qualitatively 
very similar. 

We use the following parameters: the 6-quark pole mass rrib = 4.8 GeV; the renormal- 
ization scale of the QCD tensor current u = rrib', the initial scale for the renormalization 
group evolution Hh = mi,; the renormalization scale and final scale of the renormaliza- 
tion group evolution /i = 1.5 GeV. This is also the (hard-collinear) scale at which all 
other scale-dependent quantities such as meson light-cone distribution amplitudes and 
the scale-dependent decays constants fs, fhii. are evaluated. The strong coupling is ob- 
tained from as{mi,) = 0.215 by employing 2-loop running (Ag|-^) = 323.6 MeV), which 
gives as(1.5GeV) = 0.359. The pion and p meson parameters are /^r = 130.7 MeV, 
fp\\ = 209 MeV, fp± = 150 MeV, and the second Gegenbauer moment is assumed to be 
= 0.1 for the pion and the distribution amplitudes of both, the longitudinal and 
transverse p meson. The B meson mass is = 5.28 GeV and the decay constant 

= fB/K{1.5GeV) = 200 MeV. We assume the model dHll) for the B meson distribu- 
tion amplitude and Xb = 0.35 GeV. This is somewhat smaller than the value 0.46 GeV 
suggested by QCD sum rule calculations jSl]. Allowing As to vary from 0.25 GeV to 
0.5 GeV implies that the value of A^ is the single most important uncertainty in the final 
numerical calculation. The SCETj form factors C,a{E) are defined in the physical form 
factor scheme through full QCD form factors according to (j24}. The full QCD form fac- 
tors needed for this definition are taken from the light-cone QCD sum rule calculations 
IHS] including the parameterization of their dependence. On the basis of this input we 
can compute the remaining seven pion and p meson form factors using We relate 
hadronic to partonic variables by first eliminating E through E = xmh/2 in the coeffi- 
cient functions. The energy fraction x is then interpreted as x = 1 — /m^ = 2E/mB, 
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when we plot hadronic form factors as functions of or hadronic energy E. 



5.1 Symmetry-breaking corrections to form factor ratios 

In the absence of radiative and power corrections, the factorization formula (0) im- 
plies parameter- free relations between form factors j7], since only ^a{E) appears on the 
right-hand side, which cancels in ratios. These relations receive corrections, which are 
calculable at leading power in the l/m;, expansion given the above-mentioned input pa- 
rameters |H] . The seven relations between the total of ten pion and p meson form factors 
are obtained from the two relations (j2Sl), which do not receive any perturbative correc- 
tions, and the five relations that follow from ()25|) by dividing through the appropriate 
^f^. For instance, the second and third equations of imply 

mB + mvTi{E) _ ^ ^ ^ + my j^^Bi), j^^E±{t,E) 



with C^_^\t, E), Cq^y {t, E) the combinations of coefficient functions defined in (jll9p . 
Similar relations follow for the other form factors. The second term on the right-hand 
side equals the hard spectator-scattering term piSj) divided by the appropriate C,^^{E). 
Putting together ()26|) . pi6p and ()117p we obtain the form factor ratios including the 
new next-to-leading order (and resummed) correction to the spectator-scattering term. 

The result of this computation is shown in Figure El for the various form factor ra- 
tios. The ratios are normalized such that in absence of any radiative corrections they 
equal 1 for all q^. Our final result, which includes R to order as and the spectator- 
scattering term to order as well as the summation of leading logarithms to all orders 
is shown as the solid (black) curves. To display the size of the various contributions 
to the complete result, we also show the result following from neglecting the spectator- 
scattering term (dash-dotted (red) curves), and from evaluating the spectator-scattering 
contribution in leading-order (long-dashed (blue) curves), which corresponds to the pre- 
vious results [H]. As has already been discussed in Section |3] the new NLO correction 
always enhances the symmetry-breaking effect. The correction from R{E) in ()121|) is al- 
ways smaller than the spectator-scattering contribution. In fact, it is even smaller than 
the NLO spectator-scattering term, despite the fact that the latter is formally of order 
a^. Overall, the deviations from the symmetry-limit range up to 40%, which is signif- 
icant but not anomalously large given that the typical scales involved are in the range 
of 1.5 GeV. The theoretical uncertainties in the relative NLO spectator-scattering term 
have been discussed before in Section |3J The more important unknown factor resides in 
the normalization of the tree contribution which involves the product 

XbUE) ^ ' 
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Figure 6: Corrections to the i? — > tt and B ^ p form factor ratios as a function 
of q^. The ratios equal 1 in the absence of radiative corrections. Sohd (black) line: 
full result, including NLO and resummed leading-logarithmic correction to spectator- 
scattering; Long-dashed (blue): without NLO+LL correction to spectator-scattering; 
Dash-dotted (red): without any spectator-scattering term. Dashed (black): QCD sum 
rule calculation. The lower right panel shows the two form factor ratios that equal 1 at 
leading power. For comparison, the QCD sum rule results for these two ratios are shown 
(upper (blue) line refers to Ai/V , lower (black) line to T2/T1.) 
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of hadronic parameters. We estimate the theoretical errors of the factors to be around 
15% ifB), 15% {fp±), 10% {{v-^)m), 30% (As) and 15% i^a{E)), so it is clear that the 
curves in the figure are affected by a significant normalization uncertainty. In particular, 
adopting the QCD sum-rule result = 0.46 GeV rather than 0.35 GeV decreases the 
deviations of the form factor ratios from unity by about 30%. 

It is instructive to compare this result for the form factor ratios with the QCD sum 
rule calculations. The corresponding sum rule ratios are shown as dashed (black; black 
and blue in the lower right panel) curves in Figure IHl One notices that the sum rule 
calculation satisfies the symmetry relations remarkably well - the ratios are in general 
closer to 1 than predicted on the basis of the heavy-quark limit corrected by radiative 
and spectator-scattering effects. There are also significant differences concerning the 
sign of the correction similar to those observed already in [H]. It is unclear whether 
the differences between the sum rule calculations and those based on the heavy-quark 
limit are due to 1/mb power corrections or ununderstood systematics of the sum rule 
calculation (see the discussion in jSEl)- For instance, the sum-rule result for the ratio 
involving mB/{2E)T2 — T3 has changed from about 0.7 to almost 1.2 with the update 
[53] of the form factor calculations. This may not be surprising, since the ratio involves 
cancellations between form factors and may be particularly sensitive to the uncertainties 
of sum rule calculations. In such cases the SCET calculation of form factor relations is 
probably more reliable than the QCD sum rule method. In general, the comparison of 
the two methods leads to the conclusion that the theoretical calculations of form factors 
with QCD sum rules are affected by considerable uncertainties until the systematics of 
and discrepancies with the heavy-quark limit are better understood. 



5.2 Radiative vs. semi-leptonic decay 



Factorization calculations of radiative and hadronic two-body B decays involving only 
light mesons (and leptons) make use of the form factors at maximal recoil. It is therefore 
of interest to investigate the short-distance corrections at x = 1, i.e. E = 771^/2 or 
= 0. In addition to the exact relations (|23p. the first and fourth relations of (|25|) also 
degenerate to 



FF 



(123) 



as a consequence of the equations of motion. This leaves only two interesting ratios, 
namely TZi and 7^2 defined in (jl2H) . The last ratio in (j25p involving T2 and T3 can be 
obtained from TZi replacing = /_|_ by 

At X = 1 we obtain the analytic expressions (assuming the asymptotic form 0a/ ("w) = 
6vv for the light-meson distribution amplitude) 



7^l 



1 + 7^ 

4:71 



V 



--In h - 

3 mfe 3 



8 , V 8 , 9 m;, 4 

- In \yv - 

3 mft 3 /i 3 
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1 + 0.046 (Rt) - 0.165 {l + 0.540 (NLO spec.) - 0.01 (61^)} 



0.794, 




8 u 4] 47ra,(/i) /b/m± 



3 m;, 3j 3 MbXbC 



( 



1 + 




1 - 0.023 {Rt) + 0.086 {l + 0.418 (NLO spec.) + 0.03 {61^}} 



1.102, 



(124) 



where 6^^^ denotes the small effect form the renormalization-group summation. The nu- 
merical results refer to the pion (TZi) and p meson (7^.2) with the parameters as specified 
above. The ratio TZi with = /+ replaced by and pion parameters replaced 
by p meson parameters gives 0.707 instead of 0.794. For comparison the QCD sum 
rule calculation [^2] gives TZi = 0.96 (1.02 for p and the relation involving T2, T3) and 



The factorization approach allows us to make predictions for the exclusive radiative 
B decays ^ M7 ETj and ^ Ml+l~ fSH]. The decays B ^ p'j together with 
B — s> ii'*7 are particularly interesting, because they may lead to a determination of the 
CKM matrix element \Vtd\ or constrain flavour non-universality in penguin transitions. 
The main limitation turns out to be the poor knowledge of SU(3) flavour symmetry 
breaking in the ratio of tensor form factors Tf (0)/T/^ (0) I40j . In j^H] it was therefore 
suggested to take the ratio of the B ^ p^ to the differential semi-leptonic B — > plu 
branching fraction, which avoids the problem of SU(3)-breaking, but introduces the 
ratio Ti/V of p-meson form factors at = 0. The method relies on normalizing the 
B ^ P'-f rate to the differential decay rate 



near cos^ = 1 (with 6 the angle between the neutrino momentum and the B meson 
momentum in the lu center-of-mass frame) and q"^ = 0. The angle cut has the effect 
of isolating the negative helicity form factor which has a simple expression in the 



7^2 = 0.95. 



T{B plu) 
dq'^d cos 9 



oc (1 + cos ey if! + (1 - cos 9f {Hi + 2H^) 



(125) 
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heavy-quark limit. Neglecting quadratic effects in the light meson mass, 



niB 




TUB 

hence the ratio of branching fractions involves 



V{q'), (126) 



~ ml [ ml) ^ mini' 



(127) 



Assigning a 60% uncertainty to the spectator-scattering contribution to 7^2 "we obtain 
l/7^2 = 0.82 ± 0.12 to be compared with the QCD sum rule value l/7^^ = 1.11 ± 
0.22, where we assigned a 10% uncertainty to the calculation of The disagreement 
between the two numbers is unfortunate and should be resolved. Assuming the result of 
the calculation in the factorization approach, we obtain a 10% uncertainty on iVJd/Vu^l 
from the form factor ratio in the method proposed in [32]. This does not include an 
uncertainty from power-suppressed effects. 

The tensor-to- vector ratio lZ2{q^) also appears in the forward-backward asymmetry 
in the electroweak penguin decays B —>■ Ml^l~. The complete calculation of the decay 
matrix element divides into "factorizable" and "non-factorizable" contributions [HI 123 • 
In this terminology, "factorizable" contributions are related to the heavy-to-light form 
factors, and hence are the relevant ones here. Inserting 

Ti T2 
(ms + my) — = {ms - my) — = mB'R.2 (128) 
V Ai 

into Eq. (75) of [2Hi; we obtain the differential forward-backward asymmetry [H] 



dA 



FB 



dq^ 



oc Re 



2mhMB „efr^ / 2n 



Cg + F(g") + — ^-^ Cf 7^2(g') + non-fact, terms] . (129) 



Since the dependence on is mainly through 7^.2 (g^)/^^ it follows that the increase of 
^2(9^) by several percent due to the next-to-leading order spectator-scattering correction 
implies an increase in the position of the asymmetry zero in approximately the same 
proportion. 



5.3 Hadronic decays 

The jet-function computed in this paper also appears in the next-to-leading order correc- 
tion to spectator-scattering in hadronic two-body B decays to light mesons. We outline 
this effect by the example of 5 — > vrvr decays, keeping the discussion short, since the 
NLO correction is not yet completely available. 
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The factorization formula reads PP 



{TXTimB) = /+(0) * + * fB(t>B+ * * U4>n, (130) 

where Qi denotes an operator from the effective weak Hamiltonian, and the formula 
holds up to power corrections. The second term describes spectator-scattering. Its 
short- distance coefficient is a convolution T/^ = C/^ * Jy, where C/^ is the coefficient 
function of a generalization of a B-type operator that takes into account the second 
pion. Since the pion that does not pick up the spectator anti-quark from the B meson 
decouples at the hard-scale, the physics at the hard-coUinear scale is exactly the same as 
in the 5 — > vr transition. Hence the jet-function in hadronic decays equals J\\ [21], which 
has been computed above. Note that this implies that the strong rescattering phases 
are all generated at the hard scale mi, (at leading order in the heavy-quark expansion), 
since the jet-function is real. 

Spectator-scattering is particularly important for decay amphtudes of the "colour- 
suppressed" final state 7r°7r°, because the colour-suppression is absent for the spectator- 
scattering term. The situation is opposite for the colour-allowed final state vr"?!"*". Both 
amplitudes are relevant to 7r~7r°. In the following we shall therefore focus on the coef- 
ficient a2{7in) that describes the colour-suppressed tree amplitude. We emphasize that 
a complete NLO calculation of spectator-scattering requires the calculation of the hard 
coefficient C-^ as well. The remarks below must therefore be understood as preliminary. 

Following the notation of jH] (Eqs. (35) and (47)) we write 

r \ Ci Ci as{fi)CF -.^ , . 

ao TTTT = Co -\ Vo TT 

^ ^ Nc 47r ^ 



i7*"2(7r7r)/|| +iJ*^=^(7r7r) 



(131) 



where now /i should be chosen of order mi, and is a hard-collinear scale assumed to 
be = -\/A/i/i with = 0.5 GeV. The Ci are Wilson coefficients from the effective 
weak-interaction Hamiltonian, V2{'n') is a vertex correction, and H2'"'^ {tttt) + (tttt) the 
spectator-scattering term at tree level, which we separated into a leading-power ( "tw2" ) 
and a power-suppressed "chirally enhanced" ("tw3") term. The new ingredient in this 
formula is the factor /y, which equals 1 in the absence of the NLO correction to the jet- 
function, and is given by including the correction. Exactly the same modification 
applies to the spectator-scattering contribution to ai(7r7r) and the leading-power pieces 
of the penguin amplitudes. Numerically, with parameters defined in [l^, we obtain 



[0.11 ■ 1.37]Htw2 .^ij + [0.07]Htw3 (defauh) 
[0.29-1.57] [0.17] '„3 (S4) 



a2(7r7r) = 0.17- [0.17 + 0.08i ]v2 + 

0.22(0.18) -O.OSz (default) 

0.64 (0.47) - 0.08« (S4) ^ ^ 
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Scenario 


deiauit, LU jet 


1 i? 1 J_ ATT / \ • J 

default, NLO jet 


S4, LO jet 


Ci A ATT / \ • j_ 

S4, NLO jet 


ai(7r7r) 


(J. 99 + U.U/z 


n no 1 n no ' 

0.98 + 0.02« 


n 1 n no ' 

0.88 + 0.02« 


n "1 1 n no ' 

0.81 + 0.02« 




0.18-0.08Z 


0.22 - 0.08i 


0.47 - 0.08i 


0.64 - 0.08i 


Br(i?o ^ TT+TT-) 


8.86 


8.62 


5.17 


4.58 


Br(BO ^ ttOttO) 


0.35 


0.40 


0.70 


1.13 


Br(5^ ^ tt-ttO) 


6.03 


6.28 


5.07 


5.87 



Table 1: Tree amplitude coefficients ai and 02, and the CP-averaged vrvr branching ratios 
in units of 10~^ in the default and S4 scenario of jl^ showing the effect of the NLO jet 
function correction. 

The various terms and factors correspond to those in (jl3Hl and we show the numbers 
for the default parameter set and the set S4 that provides a better overall description 
of hadronic two-body modes. Due to the near cancellation^^ of the tree term with the 
vertex correction the colour-suppressed tree amplitude comes essentially from spectator- 
scattering. The factors 1.37 and 1.57 show the effect from the NLO correction to the 
jet-function. In the final line the number in brackets gives the result from ^Tj, the 
unbracketed number corresponds to including the new jet-function term. To illustrate 
the implications of these results, we show in Table^the CP-averaged B —>■ vrvr branching 
fractions corresponding to the four cases (default vs. S4, with and without NLO jet- 
function correction). For simplicity, we only consider the NLO jet function correction 
to the tree amplitudes ai(7r7r) and a2{7i7i) (colour- allowed and colour-suppressed), but 
not to the penguin amplitudes, since this gives the dominant effect (and the results are 
preliminary anyway, see above). It is clearly seen that the NLO correction to spectator- 
scattering can have a significant effect. The enhancement of the colour-suppressed tree 
amplitude brings the theoretical computation in better agreement with data, since the 
large 71^71^ rate and the small 7T~^7r~ to 7r~7r° ratio favour a large colour-suppressed tree 
amplitude [Hj. We do not discuss the direct CP asymmetries, since we expect the still 
missing NLO hard correction to spectator-scattering (which includes a new source of 
rescattering phases) to be the more important factor. 

6 Conclusion 

Spectator-scattering plays an important role in the theory of exclusive B decays. It 
is also rather complicated, because several scales, mh (hard), y/rufji (hard-collinear), 
and A (hadronic) are involved. The development of QCD factorization and soft-coUinear 
effective theory has made it possible to formulate the calculation in terms of two separate 

^^The size of the loop correction is due to the absence of colour-suppression, which makes the tree 
amplitude small, and is therefore not an indication of failure of the perturbative expansion. 
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matching steps, in which the effects from the short-distance scales mt and A/m^A are 
calculated in perturbation theory. In previous work ^H] we began the calculation of 1- 
loop corrections to spectator-scattering effects in heavy-to-light meson form factors in the 
large-recoil region with the computation of the hard coefficient functions. In this paper 
we have completed the second step with the computation of the hard-collinear coefficient 
function, also called jet-function. Since the calculation involves the definition of various 
renormalized operators in QCD, SCETj, and SCETh, and the treatment of evanescent 
operators in dimensional regularization, we have described the technical aspects of this 
work in some length. Our results provide a check of similar results obtained by Becher 
et al. fniED]- The jet-function computed here is relevant to many different B decays, 
including radiative and hadronic B decays in the QCD factorization approach. 

The results may be summarized as follows: we find significant enhancements of 
spectator-scattering at next-to-leading order, which increase the deviation of form factor 
ratios from the asymptotic heavy quark limit, in which perturbative and power correc- 
tions are neglected. We have also included the summation of formally large logarithms 
Inmfe/A, but found this effect to be negligible compared to the full 1-loop correction. 
Despite the small scale of order 1.5 GeV involved, there is no sign that a perturbative 
treatment is not applicable. The 1-loop effects from the hard scale and the hard-collinear 
scale are about equally important, being on the order of (20-40)% (depending on pa- 
rameters), at least in the MS factorization scheme adopted throughout this work. It 
follows that the dominant theoretical uncertainties are related to hadronic input param- 
eters such as moments of light-cone distribution amplitudes and decay constants. In 
addition to the symmetry-breaking corrections to form factor ratios, we also discussed 
radiative and hadronic two-body decays. Although the jet-function constitutes only one 
aspect of the NLO correction to spectator-scattering in hadronic decays, we have seen 
that the NLO enhancement has interesting implications for final states with significant 
colour-suppressed tree amplitudes. 

We would also like to emphasize the theoretical conclusions from this calculation, 
since the form factors are up to now the only observables, for which a complete two-step 
matching in soft-coUinear effective theory has been explicitly carried out to the 1-loop 
level in a case with spectator-scattering. The factorization arguments that lead to the 
formula (Q) rely on the demonstration that the B-type SCETi currents can be matched 
to SCETn without encountering endpoint-divergent convolution integrals, which would 
violate naive SCETn factorization ^21 El- The calculations performed here and in ^H] 
provide an explicit verification of these arguments at the 1-loop level. 

Note added: We have been informed of related work by G. Kirillin, in which he computes 
the 1-loop correction to the jet-function Jy, and to the coefficient function C]-^ . 
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A Short-distance coefficients 
A.l Change of basis 

The coefficient functions of the operators defined in (jH)) to ()10|) are given in terms of 
those defined and calculated in ^Hl (denoted with subscript "old") as follows: 

scalar : 

MAO) _ ^{AQ) 
— '-^5old' 

MAO) 

=4S-%^, (133) 



vector : 

^(^0)1 _ MA0)1 MA0)2 _ MA0)3 

MAO)Z _ MA0)2 MAO)l 

MB1)1 _ ^(Bl)3 _ 1 (c)MAO)l ^{A0)3\ 



'V ~ '-'Void ^ y^^'-'Vold "T '-^Vold J ) 
MB1)2 _ MB1)2 1 fr>MAO)l MA0)2 



'V — '-'Void "T ~ \'^'~'Vold "-"Void j ) 

1 MA0)2 
MB1)3 _ MB1)1 , J- ^(Bl)4 <-^yold 

— '-^Vold "I" 2 ^Vold ; 



^(Sl)4 _ 1 ^{31)4 

- 2 '-'V'old ^ 



7(4T-4r), (134) 



tensor : 



MAO)\ _ 1 ^(AO)l ^(A0)2 _ ^(A0)3 

— 2 Told ' — ^To\d ' 

^(A0)3 _ MA0)2 MAO)\ ^(A0)4 _ ^(A0)1 ^(A0)3 ^(A0)4 

— Ogn^id — ^To\d 1 — '-Told + '-Told '-Told ) 

MBl)l _ MBl)2 1 ^(Bl)6 1 /^{A0)1 ^{A0)4\ 

Oj. — ^Told 2 Told "T ~ l^'-Told + '-^Told ) ' 



^{51)2 _ MBl)l ^{Bl)5 C'toM 
— Oj. '-Told 



(A0)1 
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'-^T — ~2 ™<i ~ X ^ "'"'-'Told ~ '-^Told J 5 

^(Bl)4 _ ^(Bl)5 1 /^(A0)1 ^(A0)2\ 
•-^T — '-Told "T ~ l^'-Told '-"Told J i 

MBl)b _ MBl)3 _ 1 /r,^(A0)2 9^(A0)3 _ ^(A0)4\ 

— '-'Told ^ l^^'-'Told "T ^'-'Told '-'Told ^ ' 

^{31)6 _ 1 ^(Bl)4 1 /^(A0)1 r,^(A0)2\ 
'-"T — 2 Told 2x V '^oW ^'-"Told J ) 

^{51)7 _ 1 /^{m)7 ^(Bl)4\ 1 /^(A0)1 MA0)2\ 

— ^ l^'-^Told ~ '-"Told J "•" 2^ l^'-'Told ~ '-'Told ) ■ K^'^^l 

In the new basis the pseudoscalar coefficients equal the scalar coefficients, and the axial- 
vector coefficients equal the vector coefficients. Furthermore x = ri-V n+p' / nib = 2E/mb 
with E the energy of the light meson. 



A. 2 Coefficients appearing in the form factors 

The five independent AO-coefficients appearing in the SCETi representation of the form 
factors (PT|) . are given by 



^(AO) 
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(136) 
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The variable E used in (1211) . (j22p is related to x through x = n_vn^p' /mf, = 2E/mh. We 
also define = as{fi), fi = fi/mh and z> = u/rrii,, where u is the renormalization scale of 
the QCD tensor current, and fi is the SCETj renormalization scale. The fi dependence 
cancels the corresponding dependence of the SCETj form factors ^a{E). The heavy quark 
mass is renormalized in the pole scheme. The five independent B-coefficients are given 
by 
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The variables E and r used in ()21|). are related to x and ^ through x = 2E/mi, 
and ^ = r. Diagrammatically corresponds to , the fractional longitudinal 

momentum carried by the transverse coUinear gluon in the B-type current operator. We 
also use ,^ = 1 — and introduced the two abbreviations 



F{y, z) = y — \ny + Li2(l — — In^ z + In^; — Li2(l — z), 

G = Li2(l - x) - Li2(l - xO - Li2(l - xO + — . (138) 



These results are obtained by taking the appropriate linear combinations of the coeffi- 
cients given in |THI. The variables (xi,X2) used there are related to (x, ^) by Xi = x^ 
and X2 = x^ (^ e [0, 1]). 



B Integrals of coefficient functions 



B.l Integration of the jet-function 

The integrals Jq dr ja{T; v,uj) of the jet-functions (jZHl), are given by 
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The analytic expressions of the integrals entering (see flH7|) ) read (setting ri/ = 4, 
T/ = 1/2 and Cj. = 4/3 and Ca = 3) 
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B.2 Convolution of Cx with the hght-cone distribution am- 
pUtude 

Because the expressions are lengthy, we only list the results for the convolution of the 
three combinations of coefficients functions in the physical form factor scheme as de- 
fined in (I25|) . and assume that the light-cone distribution amplitude are given by their 
asymptotic forms (pAii^) = Qvv. The convolution integrals read 
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with X = 2E/mh, fi = fi/mi,, = u/nih, and = as{fi)- 
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